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Abstract. We investigate the long term behavior in terms of global attrac- 
tors, as time goes to infinity, of solutions to a continuum model for biological 
aggregations in which individuals experience long-range social attraction and 
short range dispersal. We consider the aggregation equation with both degen- 
erate and non-degenerate diffusion in a bounded domain subject to various 
boundary conditions. In the degenerate case, we prove the existence of the 
global attractor and derive some optimal regularity results. Furthermore, in 
the non-degenerate case we give a complete structural characterization of the 
global attractor, and also discuss the convergence of any bounded solutions to 
steady states. In particular, under suitable assumptions on the parameters of 
the problem, we establish the convergence of the bounded solution u (t) to a 
single steady state Ut and the rate of convergence: 

|« (t) — Ut llj^pjQ) ^ {I + t)^ , ast— )-oo, 

for any p > 1, and some p = p{ut,p) g (0,1). Finally, the existence of an 
exponential attractor is also demonstrated for sufficiently smooth kernels in 
the case of non-degenerate diffusion. Our analysis extends and complements 
the analysis from [171 and many other fundamental works. 

Contents 

1. Introduction 1 

2. Weak solutions 4 

3. Optimal regularity and the global attractor IC 

4. The non-degenerate case and smooth kernels 19 

5. Convergence to steady states 31 
References 36 



1. Introduction 

In this article, we are interested in the asymptotic behavior, as time goes to 
infinity, of general aggregation models with (non) degenerate diffusion of the form 

(1.1) dtu + V ■ (ul^) ^ AA (m) , infix (0, oo) , 

where 1? = V/C * u, and $7 is a bounded domain in K'', d > 1. The kernel JC 
incorporates the sensing range and degradation for the particular population density 
u under consideration, while the term on the right-hand side of (|l.ip models the 
dispersal mechanism, such as, local repulsion. We aim to investigate (jl.ip with both 
no-flux boundary conditions and Dirichlet boundary conditions. More precisely, we 
wish to consider the situation in which the boundary dfl of the domain fi consists of 
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two disjoint open subsets Fa? 7^ and Td (possibly empty), each Ti\Ti {i e {D, N}) 
is a cr-null subset of dft and dil = Fjv U Td ; Here, a denotes the restriction to dft 
of the {d — l)-dimensional HausdorfF measure which coincides with the Lebesgue 
surface measure, if we assume that dfl is at least Lipschitz. Thus, we consider the 
following boundary conditions for (11.11) : 

^ 2) / (VA (u) - ul^) • it = 0, on Tat X (0,00) , 

' ' I w = 0, on X (0,cx)), 

and initial condition 

(1.3) W|t=o = "0 in ^■ 

The physical motivation for taking boundary conditions as in (jl.2p is clear. For 
equation (II. ip . the homogeneous boundary condition ouTd may be interpreted as 
if the species suffers extinction if say the patch C dfl where the individuals live 
is toxic. The no-flux boundary condition in ()1.2|) says that nothing can cross the 
boundary T n ■ Dirichlet boundary conditions for the population density u can also 
arise for reaction-diffusion systems in the modelling of competition between two 
population species whose interaction occurs mainly in a region where their habitats 
overlap. This gives rise to Dirichlet boundary conditions for either species on the 
whole on, or only on a part To of dil (see [44j |49j for the biological literature; cf. 
also [13l |42] , for some mathematical results) . 

We aim to give some results which allow to deduce the — L°°, and then 
the L°° — C" (ri) smoothing properties for solutions of (|l.ip - (|1.3l) assuming that 
some sort of energy estimate is apriori known in L^-norm for some finite p. The 
main tool will be an iterative argument following a well-known Alikakos-Moser 
technique combined with a suitable form of Gronwall's inequality and a refined 
ODE argument. Our estimates are much stronger than those obtained in [121 Hi] 
since our constants are uniform with respect to time and the initial data. It is 
well-known that the above smoothing properties become essential tools in attractor 
theory where they can be used to establish the existence of an absorbing set in 
C" (f2)-norm if this property can be deduced easily in L^-norm for some finite p. 
Recall that a subset V C H, where "H is a topological space endowed with a given 
metric, is called absorbing if the orbits corresponding to bounded sets B of initial 
data enter into V after a certain time (which may depend on the set B) and will stay 
there forever. If the space V is further compactly embedded in H, then the existence 
of the global attractor for a system like (|l.ip - ()1.3p follows from standard abstracts 
results (see [MlHl]). The global attractor for ()l.ip - ()1.3p encodes all the information 
about the the long-time behavior of solutions to (jl.ip - ()1.3|) departing from bounded 
sets of initial data. Our aim is to establish such a result for our system. Wc will 
consider diffusions like A{y) ~ y™, for some m > to* > 1 (the subcritical case). 
This scenario is in complete agreement with biological observation, that only in the 
case of degenerate diffusion there are solutions of ()l.ip - (|1.3p which have compact 
support, steep edges, and a constant internal population density (cf. [53]). We 
shall also emphasize the role of boundary conditions, see (jl.2p . and how they affect 
various dissipative estimates. In fact, when a (Tu) = 0, and to* = 1 + I/7 — 2/(i, 
for any 1 < 7 < d/2, we extend the results in [Ij to show the existence of a 
global attractor, bounded in C" (fj) , for some a € (0, 1) , as long as m > to*. In 
the critical case to = to*, the same result is valid provided that for J^uodx = M 
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(see (|2.ip below), we have M < Mc, where Mc is the critical mass estimated in 
[TTI Theorem 7 and Proposition 3]. However, when a (Td) > it appears that we 
cannot recover the critical exponent m* as above. We can only show the existence 
of the global attractor, bounded in C" (ft) , provided that m > = 2 in this 
case. This appears to be due to loss of the conservation of mass property in the 
associated non-degenerate problem (i.e., if A' (y) > e > 0, Vy > 0) which becomes 
crucial in obtaining uniform apriori estimates for the degenerate problem 
(|1.3|) . Finally, as in |17j we give optimal assumptions on the interaction kernel 
/C which include important cases of interest, such as, the Newtonian and Bessel 
potentials for d> 2. It is interesting to note that when K, (x) — Af (x) is Newtonian 
Af (x) — Cd |a;p^'*^^^ , for d > 3, or Af [x) — — Clog(|a;|) , for d = 2, we recover 
in (|l.ip - (jl.3p the classical parabolic-elliptic Patlak-Keller-Segel (PKS) model of 
chemotaxis: 

,^ ^, f (9tu + V • (uVc) = AA (u) , in X (0, oo) , 

^ ' ' \ c = TV * u, in X (0, oo) , 

that is, the second equation in (|1.4p reads —Ac = u. We could not find a proof 
for the existence of the global attractor, and its properties for the PKS model 
(|1.4p in the literature. Our results cover this important case as well. For the 
precise statements of the results, we refer the reader to Sections 3-4. Another basic 
system of equations modeling chemotaxis was established in the 1970's by Keller 
and Segcl [34, 35, 36 . In these systems, the population concentration u satisfies the 
first equation of p.4p . while the chemotactic agent c satisfies instead the parabolic 
equation: 

(1.5) rdtc — Ac — -I- I32U, in Q, x (0, oo) , 

for some real constants fii, /Sj and r > 0. The latter system is usually referred in the 
literature as the parabolic-parabolic Patlak-Keller-Segel model of chemotaxis. The 
mathematical and biological literature, concerning primarily with the qualitative 
properties of the solutions to the parabolic-parabolic PKS (and some of its general- 
izations) , is quite extensive and much of the work before 2003 is largely referenced 
in the survey papers by Horstmann [40l [41] . More recent results pertaining to the 
long-term behavior of solutions, in terms of global and exponential attractors, of 
the parabolic-parabolic PKS model can be found in [1] [21 |3l [25] . The issue of the 
convergence to single stationary states as time goes to infinity is also addressed in 
[221 [53] . Finally, it is worth observing that the parabolic-elliptic PKS model p.4p 
follows for an appropriate choice of the parameters Pi,/32, not only as a formal 
limit but also rigorously as r — 0+ , from the parabolic-parabolic PKS model [6] . 

Now, we wish to provide the reader with some further background on the above 
system (|l.ip - (|1.3p . The whole issue of well-posedness of weak solutions for equation 
(|l.ip with no-fiux boundary conditions on a bounded convex domain of class , and 
sufficiently smooth interaction kernels /C (see below for the precise assumptions), 
was estabhshed in [121 [17] (see also [7] for related results). When K. is not smooth 
enough, finite time blowup of some solutions can occur (see, for instance, [131 [HI 
[TSl [TBI [17]). However, it is worth emphasizing that in population dynamics, the 
non-local effects are generally modelled with smooth, fast-decaying kernels. When 
problem (|l.ip - (|1.3p is uniformly parabolic (i.e., A (y) > e > 0, for all y > 0), we can 
improve our analysis from degenerate diffusion. More precisely, we give a complete 
characterization of the global attractor As in this the union of all unstable 
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manifolds generated by all equilibria (steady-state) solutions of the non-degenerate 
aggregation equation (see Section lU Theorem 14. 8|) . At this point one could argue 
that the long-time behavior of the system (ll.ip - (jl.3|) with nondegenerate diffusion 
{A (y) > e > 0) is properly described by the global attractor. However, it is well- 
known that the global attractor can present several drawbacks, among which we 
can mention that it may only attract the trajectories at a slow rate, and that it 
may miss important transient behaviors because the global attractor consists only 
of states in a final stage. This phenomenon is already present for models of pattern 
formation in chemotaxis (see |52|). Another suitable object which contains the 
global attractor, and thus is more rich in content than the global attractor is the 
so-called exponential attractor (see [45l[52]; cf. also Section |4]). In Section |4] we 
show that the problem (|1.1[) - (I1.3I) with nondegenerate diffusion {A [y) > e > 0) 
admits also an exponential attractor (and as a result, the global attractor Ae 
is finite- dimensional), globally bounded in C" (fi), for some a G (0,1), provided 
we assume, in addition, that the kernel K. is sufficiently smooth at the origin, i.e., 

IC e W^24 (0)) . 

Here Bi (0) C M'* is the ball centered at the origin and radius equal to one. 

In the final Section [SJ we also discuss the convergence of any bounded solutions 
u (t) of the non-degenerate aggregation equation to single steady states, provided 
that $ is a real analytic function on R+, where $ (y) :— A [y] /y and $ (0) — 
$ (1) = 0. Moreover, we also establish the convergence rate of the bounded solution 
u it) to a single steady state u^, : 

\\u (t) — w*||ip(n) ^ (1 + i) , as t — > oo, 

for any p > 1, and some p = p{u^,p) G (0,1) (see Theorem I5.1|) . We refer the 
reader to Section [5] for the precise assumptions, statements and further discussions. 

As noted in |53^, the system (ll.l[) - (|1.3p exhibits interesting coarsening dynamics 
whose behavior is similar to another well-known (non-biological) model, the so- 
called nonlocal Cahn-Hilliard equation, which also exhibits behaviors in which small 
localized clumps form and merge into larger clumps over time. The latter equation 
has also been recently studied in [ini HH EH IMl EHl US] . We refer the reader for 
more related results concerning the nonlocal Cahn-Hilliard equation to [5H] , where a 
complete characterization of the long-term behavior is also given for this equation. 
Our proofs for the existence of attractors and their properties will explore various 
connections which exist between the aggregation equation (|1.1|) and the non-local 
Cahn-Hilliard equation (see [28], and references therein). 

2. Weak solutions 

We begin with some basic notations and preliminaries. Throughout the section, 
C > will denote a generic constant, while Q : M+ — )• R+ will denote a generic 
increasing function. All these quantities, unless explicitly stated, are independent of 
time, an approximation parameter e > and the initial data. Further dependencies 
of these quantities will be specified on occurrence. In the sequel, our investigation 
will be mainly divided into two cases: 

• (i) aiTu) = 0, i.e., F^i is empty, and 

• (ii) (7 (r^i) > 0, i.e., r^i is an open set of positive u-measure. 
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where we recall that a denotes the restriction to d'^ of the (d — l)-dimensional 
Hausdorff measure. 

Following [m Section 2], [T71 Definition 1] we make the following assumptions 
on A and /C: 

HI: The domain f7 C K'', d > 1, is convex and of class C . 
H2: A E C^[0,oo), A{Q) — 0, and there exists constants Ca,Ca > such 
that 

CAy"'''<A {y)<CAy"'-\ Vy > 0, 
for some ni > m,, with 

i + i-i>i, iffT(rc)-o, 
2, iffT(rD)>o, 

for any 1 < 7 < |. For a {Td) = 0, is the critical exponent defined in 
[TTI Definition 5 and Lemma 10]. 
H3: K, e W^^l {W^) n C3(R''\ {0}) satisfies the following assumptions: 

(i) K, is symmetric, K, {x) — k (|x|) and k is nonincreasing. 

(ii) k (r) and k (r) /r are monotone on r e (0, S) , for some (5 > 0. 

(iii) \D^IC{x) \ < C|xr''~\ for some C > 0. 

Note that since the function k in condition (H3) is nonincreasing, the nonlocal 
term in (|1.1[) models attraction. Moreover, these conditions imply that if /C is sin- 
gular, the singularity is restricted to the origin, so that both the Newtonian and 
Bessel potentials for d > 2 are included in our analysis. Finally, assumption (H2) 
implies that our problem ()l.l|) - (jl.3|) is subcritical in the terminology of [171 Defi- 
nition 6], and note that the total population u is preserved in time. In particular, 
there holds 

(2.1) {u (t)) ( u (t, x) ds = {uo) , for aU t > 0. 



Here, denotes the d-dimensional Lebesgue measure of fl. This property is 
lost whenever F^i is nonempty, ct(F£i) > and A is non-degenerate such that 
A' (s) > £ > 0. 

Let us recall some important properties of the kernel JC, proven in |17[ Section 
1.3], which imply a useful number of estimates in weak L^''°°-spaces, with quasi- 
norm 

\\u\\^^,^ = sup(aPAa (u))^^^, 

Q>0 

where Aq {u) := \{u > a}\ is the distribution function of u (see [T71 Lemmas 1-3]). 

Lemma 2.1. Let K, satisfy assumption (H3), (i)-(iii). 

(a) There holds: V/C € Lt^-°" (f7). Ifd>3, then K e L^^°° {n) . 

(b) Let 1? = V/C * u. Then, for any \ < p < 00, there exists a constant 
C = C (p) > such that 



,00 



(c) Let d > 3. Suppose that 7 G (l,d/2) . Then JC G ^/oc (^) */ '^'^'^ onZy if 



D^IC e L];^^ (Q). The same holds for V/C G Lf/J"' ' {ft) . 
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Remark 2.2. The above lemma shows that every admissible kernel /C, that satisfies 
(H3), (i)-(iii), is at least as integrable as the Newtonian potential. 

As in |12[ Definition 2.1] and [T71 Definition 3], our notion of bounded, nonneg- 
ative weak solutions to (|1.1|) - (|1.3[) is as follows. By (il) we denote the space of 
all functions u e {il) such that m = on Td, if cr (r^)) > 0, and by (H}j {il))* 
the dual of H]~, (ft) . Below, the space H]-) (ft) is replaced by (17) whenever 
a (Fd) = (i.e., if Tn is empty). 

Definition 2.3. Let T > be given, but otherwise arbitrary. A function u : 
17 X [0,T] -> [0, oo) is a weak solution of hl.l \) -hl.3 \) if 

(2.2) u e (17 X [0, T]) , A{u)e (O, T; B\, (17)) , 

9tu e 1^0, T; {Hi (17)) , wV/C * m G (17 x [0, T]) 

and the following identity holds: 
(2.3) 

(5tu(t) ,w)^jji^(fj))*^^i^(f^) + J VA{u (i)) ■Vw-u{t) (V/C*u(i)) -Vwdx = 0, 
for all w e Hi (17), for almost all t € [0, T] . 

Remark 2.4. Note that by u e C(0, T; (iJ|j (17))*) such that the initial 

condition in (|1.3p is understood in the weak sense of {Hi (17))*. In fact, in [HI 
Section 2] it is shown that u G C (0, T; LP (17)), for all 1 < p < oo, provided that 
< uq £ L°° (17) , so that the initial condition is also satisfied in the i^-sense. 

We have the following result concerning well-posedness of the system (ll.ll) - (|1.3p . 

Theorem 2.5. Let the assumptions (H1)-(H3) he satisfied in both cases (i) + (ii) 
(i.e., a {V d) > Oj, and assume < uq G L°° (17) . Then there exists a unique 
(global) nonnegative solution to problem i ti. in the sense of Definition[ 

Moreover, every weak solution satisfies the following dissipative inequality: 



2 



(2.4) £{u{t))+ I u{s) V<P {u{s)) -VIC*u{s) dxds < £ {uq) , 

Jq Jn 

for all t > 0, where 

(2.5) £{u{t)):= j <^{u{t))dx-]- ( [ u{x,t) JC {x - y) u{y,t) dxdy, 

$ is strictly convex on (0, oo) such that $ {y) := A (y) /y and $ (0) = $ (1) = 0. 

Proof. We only briefly mention the main steps in the proof. The uniqueness of weak 
solutions follows from [121 Theorem 2.4] (cf. also [ITl Theorem 3] for a more general 
result) with some minor modifications (see Lemma 12.71 below) since the boundary 
terms involving convolutions with the kernel /C vanish on To, if cr (Tu) > 0. 

Step 1 (Local existence). The local existence result can be carried out in a 
standard manner, by first regularizing A {y) with {y) = A{y) + ey, Ve > 0, /C 
by a sequence of smooth kernels J^K. (where is a standard moUifier), and then 
the initial data (0) = uo.e G C°° (17) n (17) , followed by passage to limit as 
£ — )• 0+ in the corresponding regularized problem P^: 

dtu^ + V • {u^VJeK, * Ue) = AAe {ue) , in 17 X (0, oo) , 

(2.6) { {VA,{u,)~u,VJ,IC'i'U,)-lt = 0, on Fat X (0,oo), 
= 0, on To X (0, oo) , 
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(see, [m Theorem 2.13], [iTl Section 3]). This is achieved by proving a series of 
uniform (in e > 0) estimates for the approximate solutions Ue (t) , in particular, by 
establishing the uniform bound in L°° (f2)-norm: 

(2.7) sup||uat)|lLoc(j,) <C, 

for some positive constant C which is independent of e > (and even the times t, T). 
This norm gives further uniform estimates (with respect to £ > 0) for solutions in 
the normed spaces of ()2.2|) . More precisely, we can obtain the following bounds: 

(2-8) II^K)|li2(o,T;//l,(^l)) 

(2.9) ll^t"llL2(0,T;(Hi,(O))*) ^ 

(2.10) \W'^l^*u\\mn.[o.T])<C, 

for some positive constant C independent of e. Moreover, for every e > it holds 

C 

(2-11) ll"e|lz,2(o,T;Hi(0)) ^ 

Thus, on account of the above uniform estimates the sequence of solutions is 
precompact in {D, x [0, T]) , for all 1 < p < oo, (cf. [17l Lemma 9], [H]), and we 
can pass to the limit as e — > in the nonlinear terms in ()2.6|) . In particular, it can 
be shown that 

A {u, {t)) A{u (t)) weakly in (O, T; H^, {Q)) 

and 

(2.12) MeV Je/C * Ue -> uV/C * u weakly in {VL x [0, T]) 

(see, [171 Theorem 1]). Let us now comment how to get the estimates (|2.7|) - (|2.1ip 
in each of the following cases: a {To) = or cr (Fd) > 0. 

• Case (i): a (Fd) — (i.e., F^i is empty). The argument for deducing (|2.7p 
relies on showing that the L^-norm of [t) (see (|2.ip ) apriori controls the 
LP-norm for 1 < p < cx), and then that the latter norm controls the L°°- 
norm of (t), see (|2.7|) (cf. also [17, Lemma 8]). The Lemma [2^1 is crucial 
for the proof of the uniform bound ()2.7|) . Under the assumptions (H1)-(H3) 
the bounds (|2Jl) - (|2lT|) were obtained in [ITl Section 3] (cf. also [12 ). 

• Case (ii): a {Yd) > 0. The arguments from [T71 HI] leading to estimates 
(|2.7p - (|2.7l) on intervals (0, T) , for some T = T (p) > 0, seem to break down 
if T]j is nonempty and a (To) > since there is no conservation of mass 
(|2.ip for (|2.6|) in this case; hence, the L^-norm of (t) is not controlled 
apriori. However, even in this case, the same proof of obtaining local-in- 
time estimates for (t) in [171 Lemma 8] applies, but we need to use the 
following inequality 

(2.13) \\u, mUn) < 2^ (IK". W - + l^l) ' fo'- any k > 0, 

instead of the usual one [13 Ineq. (25)]. Then, by arguing as in [T71 
Ineq. (27) and Lemma 8, Step 2] (note that boundary terms involving 
convolutions with the kernel K, vanish on Tu, if a (F^i) > 0; also now the 
last constant on the right-hand side of [El Ineq. (27)] depends only on 
k,p, and not on L^-norm of (0)), we also obtain a uniform in e > 0, 
i°°-bound for (t) , t & [0, T] , for some T = T(p)>0. 
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By virtue of these observations, the proof of the local existence argument in both 
cases (i)+(ii) goes exactly as in [l71 Section 3] (cf. also [12l Section 2.2]). A proof 
of the energy inequality (|2.4I) can be found in 12, Lemma 5.1], JTT, Proposition 1]. 
Continuation of the weak solutions is a straightforward consequence of the local 
existence theory (in Step 1) and the proof of [IJl Theorem 4]. 

Step 2 (Global existence). 

• First, let a (Td) = 0. Since the problem is subcritical by assumption (H2), 
the global existence result follows from |T71 Lemma 10 and Remark 9] which 
shows that (|2.7p is also satisfied uniformly with respect to t,T. 

• When (J {Yd) > 0, it suffices to establish a uniform (in e > and time 
t,T > 0) L^'-estimate for (t) and to exploit an argument similar to [T71 
Lemma 8, Step 2] to deduce the uniform (with respect to t,T,e) bound 
(|2.7p . However, the key difference with respect to the case ^(r^)) = is 
that mass is not conserved for (|2.6p : hence, the L^-norm of (t) is not 
controlled apriori and we need a different argument to deduce the uniform 
bound in L^-norm for Ug (t) , see Proposition 13.11 and Remark 13.21 below. 
The latter bound requires that the critical exponent m* = 2 in (H2). 

The proof of Theorem 12.51 is now complete. □ 

Remark 2.6. In order to establish the energy identity in (|2.4|) . it would suffice to 
show that $ (u) € (O, T; H}^ (ft)) so we can test equation (j2.3|) with $ (u) and 
/C * w, respectively. At the moment, this seems unreachable for our problem with a 
degenerate diffusion function A (u) ^ m™, like in (H2). However, note that $ (u) G 
(0, T; H}, (n)) , $ e (0, T; (n)) since $ (u) - A (u), by (H2). When the 
diffusion A in ()1.1|) is non-degenerate (say, like in problem (|2.6|) ). we can establish 
the energy identity (see the proof of Theorem 12. 5p . Indeed, in this case we can 
easily check that $' (u^) e (O, T; Hj^ (fl)) , Ve > 0, on account of (1^ and ^^JTl . 
Thus, the key multiplication of the corresponding weak formulation associated with 
([2:6)) (see with $' (u,) and /C^u^ e (O, T; Hjy (fl)) , respectively, is allowed. 

Exploiting, for instance, [12t Lemma 2.6], and owing to the convexity of $e we get 
equality in (|2.4|) for the energy £^ {u^ (t)) , associated with (j2.6p . which is defined 

by 

£e{ue{t)):^ / {ue (t)) dx - I- f f {x,t) IC {x ~ y) {y,t) dxdy 

($e is the same function as above, but with A replaced by A^). 

The next lemma gives the (Holder) continuity of solutions with respect to the 



Lemma 2.7. Let the hypotheses of Theorem \2.5\ be satisfied. Let ui (t) and U2 (t) 
be any two weak solutions of jj. j|) - l77S)) corresponding to any two initial data uqi 
and Mo2, respectively. If a {T d) = 0, we further take (woi) — Mi and (1*02) — M2, 
for some Mi,M2 > 0. For all t G [0,r], the following estimates hold: 




(2.14) 




and 




-ct 
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if a (Yd) = 0. 

Proof. We briefly explain how to get p.lSp : the estimate (|2.14|) is similar. Let 
u (t) := ui (t) - U2 (t) and observe that ipi^ yields {u (t)) = Mi ~ M2 =: M12, for 
all t > 0, but M12 ^ 0, in general. As in [121 Theorem 2.4], JTj Section 2], consider 
the Neumann problem 



(2.16) 



(t) = u (<) - 7\fi2 in f7, 
V0(t) • 7^ = on r, 



for which {</) (t)) = 0. Notice that since {u (t)) = M12, ((216)) has a solution. Con- 
sider the operator An = — Aat, with domain D (An) = {ip G (fl) : (Vt/? • ~rt)^j. — 
0}. Clearly, (j) (0) = (u (0) — M12) and recall that, due to a Poincare inequahty, 

II ii2 



A-^^l^ {u-{u)) 



Next, we see that that dtcj) {t) also satisfies (in the generalized sense) the problem 
(2.17) 



Adtcj) (t) = dtu (t) in Q, 
\7dt<t) (t) • it = on r. 



and 



<j)eL°°{flx [0, T]) n C (0, T; H}^ (f7)) , V(/) G L°° (f7 x [0, T]) . 
Thus, arguing in a standard way as in |12[ (6)-(12)], we obtain for 

V{t) =l|V0(t)||i.(f2) + (Mi2)', 

the following inequality: 
1 d 



(2.18) 

where 
(2.19) 



2dt 



i^{t)^{V^{t),dt^cf>{t)) 



{dtu{t),c^{t))=h+h + h, 



h ■■= f VA {ui (t)) - VA {u2 (t)) ■ {t) dx, 
Jn 

/2 := - / [u {t)) V/C * ui (t) • (t) dx, 
Jn 

h -=- I U2 (t) V/C * {u {t)) ■ V0 [t) dx. 

< C , since Ui is bounded) , 



Since A is increasing and bounded (i.e., A {u 
from (|2.16|) we have 



L°°(n) 



(2.20) h^- {A (ui) - A {u2)) (mi - ^2) da; + M12 / ^ (ui) - A (^2) da; 
Jsi Jn 

< C (Ml - M2f . 

To bound the I2, 13 integral terms, we integrate by parts and proceed as in [1^117) . 
We deduce 



(2.21) 



/ dijcjidjK, * uidi(l)dx — "S^ / di^idjjK. * uid, 
, , Jn , , Jn 



4>dx 
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- J [{djIC * ui) nj]\di(l)f da. 



The last term on the right-hand side of p.2ip is nonpositive since ft is convex and K. 
is radially decreasing (i.e., as in |T5], we have (V/C * ui)-lt < on Tpf). Integration 
by parts in the first term gives the bound: 

/ dij(f)djlC * uidi(f)dx < ~- I (A/C * ui) I V(/)|^ dx, 
which together with (|2.2ip entails 

(2.22) l2<C [ \D^JC*ui\\\/(l)\^dx, 



for some positive constant C, which only depends implicitly on the uniformly con- 
trolled L^-norms of ui and U2- Arguing as in JL?,, (17)], by using Lemma [2.11 (b) 
and since (M12) < rj (t) , we deduce for any p > 2, that 

(2.23) I2 < Cp-n {tf-^'" and h<C Wu^W^^^^^ v (t) ■ 
Consequently, we obtain the differential inequality 

(2.24) j^Tj (t) < Cpv {tf-"'' + (<) , Vt e [0, T] . 

As in [T51 Theorem 3.3], the idea is to fix p — p (r] {t)) in an optimal way. To this 
end, choose a sufficiently large constant C > such that p = log(C/77(i)) > 2, 
for all t £ [0,T] . First, recall that 77 (t) is bounded on t G [0,r] and that the first 
term on the right-hand side of ()2.24p dominates the second one since we need the 
estimate for rj (t) small only. Thus, the second term on the right-hand side is not 
essential and we can derive the following differential inequality 

(2.25) (t) < Ctj (t) log (^-^^ , G [0, T] . 

Here, we have used the elementary inequality rj (t) < C, for the p-chosen above. 
Integrating (|2.25p with respect to i G {S, s), we obtain 

_ C . 

Passing to the limit as J — > in (|2.26p . recalling that rj {t) is continuous, we get the 
desired estimate (I2.15p . □ 



(2.26) ?7(s)<C 



Vs G (0, T] 



3. Optimal regularity and the global attractor 

In this section, we derive several uniform estimates for the solutions of the 
problem (|l.ip - (ll.2p which are necessary for the study of the asymptotic behavior 
as time goes to infinity. In a first step, we obtain dissipative estimates for solutions 
in the spaces L^, L°° and C", a > 0, uniformly with respect to time and the initial 
data. Incidently, the estimates derived below allow one also to obtain optimal 
regularity results for the weak solutions on J7 x [t,oo), for every r > 0, associated 
with the system (|l.l|) - (jl.3p . Finally, the apriori estimates will be deduced by a 
formal argument, which can be justified rigourously by means of the approximation 
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procedure devised in [121 Section 2] , [ITl by means of (|2.6p . Regardless of the type 
of approximation procedure being used the regularity properties 

(3.1) Ue £ L°° iO,T;L"° {n))r\L^ {0,T;H}:,{n)) , Ve,r>0 

are essential in order to rigorously perform these computations. To this end, we 
shall only perform our (formal) computations to the original system (|l.ip - (|1.3p . for 
the sake of simplicity. 

The (uniform) dissipative L^'-estimate when a{TD) > is different than the 
estimate when a (Tn) = (in this case, it was obtained in [T71 Lemma 10 and 
Remark 9] ) , since in the former case there is no conservation of mass in (|2.6p . It is 
given by the following 

Proposition 3.1. Let the assumptions (H1)-(H3) he satisfied and assume that 
(7{T£)) > (Case (ii)). Then there exists constants C^:,C+ > 0, and /i > 1, 
independent of time and the initial data, such that every weak solution of 
lll.3\) satisfies 

(3.2) ||ii(i)||^p(o) < (^^y' +[C+{^i-l)t]-T^ , forallt>T>0, 

provided that p > m + 2d/ (d — 1). The constants C*, C+ and /i can be computed 
explicitly in terms of the physical parameters of the problem. 

Proof. We first begin by noting that V/C G (ft) , for any g > d/ (d - 1) , by 
Lemma [2TT| (a). This yields 

(3.3) V/C1b,(o) e (M'^) and V/C1r.\b^(o) G {R") , 

for any q > d/ {d — 1). Next, let (p £ (il) n (57) and A denote the principal 
eigenfunction and eigenvalue of 

(3.4) - Aip = Xif in n, 
with the boundary condition 

(3.5) \/lp - it = ouTn, 
such that 

(3.6) V(^-7i +</3 = 0onrD. 

Clearly, A > and by the maximum principle (see, e.g., 23;) > in O since 
cr (Td) > 0. W.l.o.g., we may assume that = ^■ 

Set now a (y) A (y) . Testing equation (12.31) by pu^^^cp, p > 1, we obtain 

(3.7) T" / (t) ^dx ^~p [ a{u (t)) Vu (t) ■ V (w""^ (t) ip) dx 

+ p [ u{t) (V/C * u {t)) V {uP-^ [t) ip) dx 

where 

11 := —p {p — 1) a (u) vP^"^ I Vm|^ ip (x) dx, 

12 '.= —p a (u) uP~^'Vu ■ Vip {x) dx, 

13 := p Ip - 1) /n w (V/C * w) • Vu {uP'^ip (x)) dx, 

14 :— p J^u (V/C * u) ■ V(/3 (x) uP~^dx. 
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First, from (H2) it is easy to see that 

(3.8) h < -Cp (p-l) I |yy|2 ^ ^ Q_ 



n 



On the other hand, setting a (u) — a (y) ^dy > ^, and using the 

definition of (p from p.4p - p.6p . we have 



(3.9) I2 = —p / Va (u) • V(y9 (x) dx — p a (u) Aip (x) dx 

Jn Jn 

—p / a{u)'S7ip -itda — p / a (u) 'S/ ip ■ it da 



-pX / a (u) (f {x) dx 



n 



< -ex / u'^+P-^ifidx. 



since a{0) — {u ~ 0, a.e. on x (0, 00)) and (p satisfies (13. 5p . Moreover, we can 
estimate the integrals /i, /2 using Holder and Young inequalities as follows: 

(3.10) 

h^p{p-'^) (u" "^^y<^(VA^ * u)^ • (^y/^u" 2^'"v-u^ 

<p{p-l)(^j wP-^+'^\'^u\^ ip{x)dx^ ' (^J^\VIC'*'u\'^uP-"'+^ip{x)dx^ ^ 

<VPip- 1) A + / I V/C * uf-^+V (2^) 
Jo 

for every 77 > 0. We must once again absorb the last term on the right-hand side of 
(|3.10p into Ii. For any q > d/ {d- 1) and s > 1 such that l/s+l/q= 1, Holder's 
inequality yields 

(3.11) Cr^ J |V/C*u|^MP"™+id2; < CJ|V/C*u||^2,(n) (^y" u^P'"'+^'>'dx^ ' 

< Crj (/C) ||M||i2,(s^) ll'"lli(p-™ + l)s(n) , 

since, from p.3p . there holds 

(3.12) ||V/C *M||^2,(f2) < ||V/C1b^(o)||^i ||M||i2,(n) + ||VA:iRd\B^(o)||^2j|M||^i(f2) 

Thus, choosing g, s > 1 in (jS.lip in an optimal way such that 2q — {p — m + I) s > 
2d/ [d— 1) , we further obtain in p.lOp by virtue of p. lip and Young's inequality, 
that 

(3.13) /3 < IIP {p -i)h + a, (/c, f7) ikiii;":+'3(n) 
<rip{p-i)h+-i I u'^-^+Pipdx + a,.ri (/c, n,p, TO, if) , 



for every 77 > 0, since p + to— 1 >p — ?7i + 3 (recall that to > 2). It follows by 
choosing sufficiently small 77 < C/2 in (|3?T0| . and 77 < (C/2)A in ((Xn)) that the 



MODEL FOR BIOLOGICAL AGGREGATIONS 



13 



integral term I3 can be completely absorbed into Ii. The term I4 can be bounded 
exactly the same way. We have 

h^p I u (V/C *u)-Vip {x) uP-^dx < e / V (x) dx + Ce (IC, n,p, m, (p) , 



for every e > and some positive constant which depends on G (fj) 
Summing up, from p.7p we deduce 



(3.14) ^ / uP{t)ipdx + CX I uP+'''-^{t)<f{x)dx<C{IC,n,p,m,ip). 



dt 

Finally, set /i := (m — \ +p) /p > 1 and 

Z {t) / uP (t) (pdx. 
Jn 

By Jensen's inequality, p.l4p yields the following inequality: 

(3.15) ^Zit) + C+{\)iZit)r <C,. 

dt 

We can now use the Gronwall's inequality (see, e.g., |5T1 Chapter III, Lemma 5.1]), 
applied to p. 151) to deduce the desired claim in p.2p . The proof of the proposition 
is complete. □ 

Remark 3.2. In the case when uq g (f2), Z (0) = limt_j.o+ Z (t) is finite, so a 
similar argument to j511 Chapter III, Lemma 5.1] gives 

(3.16) Z(i)<max|holLp(o),(^)"|, > 0. 

It is worth emphasizing again that the left-hand side of the inequality in (112) (i.e., 
CAy"^~^ < A (y), Vy > 0) is enough to establish the above assertion. Moreover, 
the above estimate p.l6p directly implies global well-posedness in the subcritical 
case for problem (|l.ip - (|1.3p for as long as cr(r£)) > (Case (u)). Unfortunately, 
we are not able to argue as in the proof of [T71 Lemma 10] to obtain the desired 
L^-estimate (similar to the case when cr(r£i) = 0, see p.43p below) since we do 
not know how to get apriori control over the L^-norm of u (t) for problem (|2.6I) . 
We emphasize again that when a (To) = 0, mass is conserved in both (j2.6l) and 
(IIIlD-([I3D so that \\u (Oll^iffi) < M \n\, for any M > such that (uq) < M. 

Next, we establish a crucial result which allows to deduce a dissipative L°°- 
estimate uniform with respect to the initial data, and which is necessary for the 
attractor theory. For this result, we consider both cases a {To) > 0. 

Lemma 3.3. Assume that there exists e > and > such that A (y) > e for 
O'll y ^ Uf Let V he a solution of the following degenerate problem: 

(3.17) dtv (i) + V • {v (t) (t)) = {v (t)) , in Q x (0, od) , 
such that V satisfies il.3\) and 

(3.18) 
If 



{S/A{v) -vV^) -it = 0, on Fw X (0, 00) , 
V — 0, on To X (0, 00) . 



< C, for all t > 0, 

(L-(O))'' 
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then there exists ji ^ {t — t) > 0, and a positive constant C = C (C, /i, e, j/e) such 
that 

(3.19) sup |k(i)|Loo(o) <Cg(l+ sup ||i;(t)|Li(f,)), 

t>r'>0 t>T>0 

/or a/Z T > r > 0. The constant C — C (/i) ^ (I > 0) is bounded if fi is bounded 
away from zero. 

Remark 3.4. Lemma 13.31 implies the boundedness of the function v {t) provided 
that (t) is a bounded vector on [t,oo). It establishes the L^-L°° smoothing 
property for the solutions of p.l7p - p.l8p . This result is analogous to the result 
obtained by Kowalczyk [37 , but it is more sharp since the function and constant 
on the right-hand side of p.lQp do not depend on the L°°-norm of the initial data 
(see, |37[ Lemma 4.1]). This is very useful if we want to produce uniform estimates 
with respect to time and the initial data. 

Proof. Step 1 (The local relation). First, we recall the following estimate which 
can be obtained exactly as in j371 Lemma 4.1, (5.1)-(5.4)]. Indeed, setting Vm — 
(\v\ — /)^, for any I > y^, multiplying equation p.l7p by the p-th power ofvi, p > 1, 
and integrating by parts using ()1.2p . it follows, after a suitable rescaling of the time 
variable t — et, that 
(3.20) 

^ / TJf+i (t) dx + I Vv^ (t) \x<Cp{p+l)( I iJf+i (i) dx + 1 

dt Jq_ P + 1 Jo \Jn 

for some constant C > independent of p, but which depends on I and e, and where 
vi (t, x) — vi (t/e, x) ,t £ (0, et) . In fact, C ^ 1/e as e -> 0. 
Secondly, set pfc = 2^^ — 1, fc > 0, and define 

(3.21) y, (t) / Uf+i (t) dx, 

for all fc > 0. Let i, /z be two positive constants such that t — fi/pk > 0, and whose 
values will be chosen later. We claim that there holds 

(3.22) yk (t) < Nk {I y) := C (m, e) {p^V ( _sup yk-i (s) + 1)^, VA: > 1, 

where C, 7 are positive constants independent of k. The constant C ~ C {ji, e) is 
bounded if /i is bounded away from zero. 

We will now prove p. 221) when 2 < d. The case d < 2 requires only minor 
modifications by using a suitable Sobolev embedding. The argument we follow is 
similar, for instance, to [28] (cf. also 27"). For each fc > 0, we define 

^ d{pk + l)-{d-2)il+pk) 
d(pfc + l)-(d-2)(l+pfc_i)' 

We aim to estimate the term on the right-hand side of p.20p in terms of the 
(J7)-norm of U,„. First, Holder and Sobolev inequalities (with the equivalent 
norm of Sobolev spaces in {VL) C LP= (17), p^ = 2d/ [d - 2)) yield 

(3.23) / tJf+'da; <( ( (zJ/)^^^ ^x) ' ( [ {vi)'-^"'-' dx^ ' 
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< c 



f 


(pfc+i) 


2 

dx + 


^ {Vlf^^" d^ 


In 









1+pfe- 



^ da; 



with Sfc := Sfcd/ (rf — 2) e (0, 1). Applying Young's inequahty on the right-hand 
side of p.23p . we get 

2 



(3.24) 



v}+P^dx < - 



v^'^-'+'dx 



for some positive constant independent of pk , and where 

rk 



^k := 



> 1 



1 - Sfc 

(in fact, straightforward computations give 9k = 2 since pk = 2'^ — 1, for ah fc). 
Note that Q^_^ (s) ^ s'^^ a.s s oo. Inserting estimate p.24p on the right-hand side 
of p.20p . we obtain the following inequality: 

2 



(3.25) 



d 



yk (t) + 





pfc+i 


-1 




4 70 





dx<c{pkV {yk-i(t) + iy 



where C, 7 are positive constants independent of fc. 

We are now ready to prove (|3.22l) using p.25p . To this end, let C (s) be a 
positive function Q : M+ — > [0, 1] such that Q (s) = for s G [O, t — ^J./pk\ , C (s) = 1 if 
s e p, -1-00) and \dC/ds\ < Pk/fJ^, if s S {t- fi/pk,t). We define Zk (s) = C (s)3^fc (s) 
and notice that 

^Zk{s)<P^yk{s) + c{s) ^yk{s). 

ds fi ds 

Combining this estimate with p.25p . p.24p and noticing that Zk < yk, we deduce 
the following estimate for Zk'. 

(3.26) -^Zk is) + C ifi) pkZk (s) < Nk {I ^i) , for aU s e [t - ^i/pk, +00) , 

for some positive constant C = C (/i) independent of fc. Integrating p.26p with 
respect to s from t — fi/pk to t and taking into account the fact that Zk (t — i^/pk) = 
0, we obtain that 

yk (t) = Zk (t) < Nk {I ^i) (1 - e-^^) , 

which proves the claim p.22p . 

Step 2 (The iteration procedure). Let now r > r > be given with t > such 
that 

(3.27) sup {yo{-'^) + l)< sup (\\v{s)\\^,^^^ + l + l\n\)<C^{l), 

S>ti=T S>ti=T ^ ' 

and define fJ. — {t —T),tQ = T and tk = tk-i — f-fPk, k > 1. Using p.25p . we have 

(3.28) _sup yk {t) < C {pkV (sup 3^fe_i (s) + 1)2, fc > 1. 

We can iterate in p.28p with respect to fc > 1 and obtain that 



(3.29) 



sup yk(^)<{Cpl){Cpl_,f{Cpl_ 



2k 
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where 

(3.30) Afc :=l + 2 + 22 + ... + 2^ 

(3.31) Bk :=k + 2{k-l) + 2^ (k - 2) + ... + 2^. 
We can easily show that {Afe} and {-Bfe} satisfy 

(3.32) Afc < C (2'^') and Bfc < C (2*^) , 

for some positive constant C independent of fc, j and /i (see, e.g., [28]). We can 
take the 1 +pfe = 2'''-root on both sides of p.29p and let k — > +oo. We deduce 

(3.33) sup ||i^/(i)|Loo(o) =^ 1™ sup(:Vfe(t))'/^'+^'=^ <C(M,e,ye)(C4), 

t>To=r' '=^+°°t>to 

for some positive constant C independent of Z, fc, iJ; and initial data. Next, we 
notice that, for any p > 1, in view of ^a + h < ^/a + -^/b, there holds 

thus, as p — )■ oo, we have from p.33p . 

(3.34) sup ||tJ (t) < 2 (C (m, e, 2/.) (C^) + /) , 

t>T' 

for any I > Ue. Rescaling back the time variables {t — te, t = t e, t = re) into 
(13.341) and p. 271) and taking I — y^, we easily obtain the desired inequality p.l9p . 
The proof is finished. □ 

We can now show the following. 

Theorem 3.5. Let < uo £ L°° (51) and assume (H1)-(H3) hold for both cases 
(i)+(ii). If a{T£,) = 0, let M > be given such that (uq) < M. Every weak 
solution of is globally Holder continuous in the cylinders [t,t ] x f2, for 

all T > r > 0. In particular, the following estimate holds: 

(3.35) sup(p(M)||^2([t^t+i],,fi) + ||u||c„/2,„([t^t+i]xo)) <C', «/CT(ri))>0. 

Finally, for every bounded subset B = B (M) C {^l), th ere exists a time t+ = 
i+ (B) > .such that 

(3.36) sup {\\A{u)\\^2(^[t,t+i]xn) + IMc/^.-ilut+iun)) ^ C", ifa{rD)=0, 

for some a > and some constant C > independent of the initial data and time. 

Proof. It suffices to prove the estimates p.35l) - p.36l) . 

Case 1. Consider first a{T£)) > 0. In this case, we can take p > d and 
p > m + 2d/ {d — 1) as large as we want in Proposition 13 . II so that 

(3.37) sup||V/C*ii(t)||^p(j-^) < C 
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since JC £ L'^/('i-i).oo (jg^d-^^ rpj^^g^ exploiting first Lemma[0(b), 

we have for \^ — V/C * u, 



sup 



(3.38) 

By Money's inequality (see, e.g 



<C(p)suph(i)|liP(j^)<C. 



(3.39) 



sup 



< C. 



L°=(0) 

Invoking now the crucial Lemma |3.3[ we deduce on account of 

(3.40) suph(i)|lioo(f,)<C. 

t>i ^ ' 



that 



We can now test equation (|2.3I) with A (u (t)), then integrate over (t, t + 1). Setting 
A (u) := J(j" A (s) ds, and exploiting p.40p we obtain for every t > 1, 

rt+l 



{A{u{t + l)) - A{u (t))) dx 



lVA(u(s))||l.(^) ds 



t+i 



\/A {u (s)) • (mV/C * u) dxds 



< 



t+i 



iivA(w(s))ir^.(^) ds+c. 

Thus, we easily deduce the first part of (I3.35p . For every 1 < p < c», we now have 
(3.41) A/C * u G LP (r2 X [l,oo)) and = V/C * u e L°° (f7 X [1/2, oo)) 

since u is bounded. Global Holder continuity results for (|l.ip - (|1.2p have been 
proven for more general classes of degenerate quasilinear equations in divergence 
form [211 [221 [321 [331 Hz]. Due to (jXiT]) . the second estimate in ([335]) (respectively, 
p.36p ) is a simple corollary of these results when m > m*, see, e.g., [5D] (cf. also 
[ia[2ll|47]). 

Case 2. The case a (Td) = is similar. Following JX Lemma 10] and Lemma 
12.11 (c), we see that for any m G [1, m*] such that 

liminfA' iz)z^-"^ > 



(3.42) 



(u-k) e L^n), the following 



and q = {2 — m) / (2 — to*) > 1, with Uk 
estimate holds (cf. 17, pg. 1703]): 

(3.43) ^ \\uk (t)||i,(^) < ~v \\uk mUn) + C {m, k,p, \\uk mimn} 

for all i > 0. for some positive constant j] > 0. It remains to note that by assumption 
(H2), condition (I3.42p is already satisfied with to = m*, so we can choose g = 1 in 
(|3.43p . Gronwall's inequality (see, e.g., [51]) applied to (|3.43p yields, on account of 
([2l3l) . that 



(3.44) 



l!«(t)|li.(o) < ||K-fc)+||^,(^) e-^' + C{M,Q,k,p), 



since \\u (t)||^i(s7') ~ \^\ {u (t)) < M (mass is conserved). Hence, for every subset 
B — B (M) C (fl), and setting — C {M, ft, k,p) + 1, we can easily find a time 
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t+ = t+ (Af ) > such that ^^A^ imphes 

(3.45) sup ||u(t)||^,(j,) <i?,. 

After establishing the bound p.44p . the same estimates p.37p - p.39p hold provided 
that V/C e L'^'°° (ri) , for q > d/ {d — 1) is satisfied. Thus, we can again reach the 
uniform L°°-estimate p.40p now on the time intervals [t+,oo). This is enough to 
get (|3.36l) once more by applying the results in [20l |47] (see Case 1 above) and to 
conclude the proof of the theorem. □ 



Define, for some given A/ > 0, 
Zdn ■- 



{uo e (n) : uo > 0, (uo) < M} , if a (Td) - 0, 
{uo e (fl) : uo > 0} ifcr(ri5)>0. 

In the previous section we have proved that system (jl.ip - (jl.3p generates a contin- 
uous semigroup S{t) on the phase space Zdn, endowed with the metric topology 
of {Hjy {ft))* , via 

(3.46) S it) : [Zdn], ^ [ZdnI , S (t) uo = u (t) , t > 0, 

where u{t) is a unique weak solution of ()l.ip - (|1.3p (see Theorem [23]) . Here [^^dat]^ 
denotes the closure of Zdn in the metric of (^H^ i^))* ■ We devote next our attention 
to the study of the long-time behavior of trajectories of the semigroup in terms of 
global attractors. We need to recall the following definition. 

Definition 3.6. We say that A C Zdn is the global attractor for the dynamical 

system {S (t) , Zdn) if the following conditions are satisfied: 
(i) the set A is a compact subset of the phase space Zdn i 
(a) it is strictly invariant, that is, S{t)A = A, for all t > 0; 
(Hi) for every bounded subset B C Zdn , 

(3.47) dist^ZoN], (t) B, A) 0, as t oo, 

that is, A attracts the images of all bounded subsets of Zdn cis time goes to infinity. 
Here, 

distizr^r,], i^^Y) = sup mf,||a;-y||[2„^]^ . 

The first main result of this section states the existence of such an attractor for 
problem (fTTIl - lfOl) . 



Theorem 3.7. Let the assumptions of Theorem \ 3.5\ hold. Then the dynamical 
system {S {t) , Zdn) associated with problem m.l\) - nTS\) possesses a global attractor 
Adn in the phase-space Zdn, which is globally bounded in C" (Tl) , a S (0, 1) and 
has the following structure: 

(3.48) Adn = ^\t=07 

where S is the set of all bounded solutions of il.l]) - nT3\) . defined for all t e R, such 
that 

(3-49) Pllc./2,.(MxTT) < 

for some positive constant C. 
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Proof. As usual, we must check that S (t) possesses a (pre)compact absorbing set 
in Zdn, and that it is closed. The first assertion is an immediate corollary of 
Theorem 13.51 Indeed, for every bounded subset B C Zdn, there exists a time 
t+ = t+ (B) > 0, such that S{t)B C X", for aU t>t+, where 

(3.50) X":^{ueC"^{n):\\u\\^.(j^^<c} 

(C is the same constant as in p.35p - (|3.36p ). Moreover, owing to Lemma [2?7l we 
also have that S (t) is a closed semigroup in the sense of [46]. Thus, due to the 
abstract theorem on the attractor's existence [46 [ Corollary 6], this semigroup pos- 
sesses a global attractor Adn, bounded in C" (fi). The characterization p.48p and 
estimate p.49p follow from Theorem 13.51 and the same abstract results. Lemma [37fl 
is proved. □ 

Remark 3.8. Of course, the global attractor Adn in Theorem 13.71 depends on 
M > since the constant C > in p.50p does. Furthermore, the entering time 
i+ > is even independent of the set B C Zdn, when a {Td) > 0, while it only 
depends on M > whenever a (Fd) = 0. 

Finally, we can extend the above result in the critical case m = and a (Td) = 

0. 

Theorem 3.9. Let assumptions (HI), (H3) and (H2) with m = m, he satisfied, 
and consider the case when T ^ is empty. Let Mc > be the critical mass estimated 
in [171 Proposition 2 and 3]. Then for every M < such that (uq) < M, the 
dynamical system {S (t) ,Zun) possesses a global attractor Am^ C Zdn, bounded 
m C" (I7) , a e (0, 1) . 

Proof. As in [17, for every j^u^dx — M < Mc, problem (jl.ip - (ll.3p is globally 
well-posed (see, in particular, Theorems 7 and 9]). The assertion in the theorem 
follows by arguing exactly as in the proof of Theorem 13.51 (Case 2) since (I3.42[) is 
still satisfied when m = . □ 

4. The non-degenerate case and smooth kernels 

In this section, we are interested in the model proposed in [7] which also takes 
into account stochastic fluctuations based on a finite number of individuals subject 
to long range attraction and short range repulsion. In this case, the density u (t) 
satisfies 

(4.1) dtu = div (VA (u) - u (V/C * u)) + eAu, in n x (0, oo) . 

The additional parameter s > models classical Brownian random dispersal in 
equation (|l.ip . which can be seen as the limit of (14.11) as e 0, see [7]. For the 
sake of simplicity, we will only consider the no-flux boundary condition for (j4.1l) . 

(4.2) (V {A (u) + eu)-u (V/C * u)) • ^ = 0, on F x (0, oo) , 
and, as before, the initial condition 

(4.3) uit=o = "0 in ^■ 

Thus, everywhere in this section, the boundary F^, where u = 0, is assumed to be 
empty (hence, F = dfl). 

As in Section 3, we have the following result whose proof is straightforward (see 

mm)- 
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Theorem 4.1. Let (HI), (H3) he satisfied and assume 
(4.4) liminfA' {z) z^-^' > 0, 

where m, > 1 is the same critical exponent used in Section 3. Ifm^ — 1, we further 
assume that each AI > satisfies M < where (uq) < M. Then there exists 
a unique (global) nonnegative solution u(t) = (t) to problem |^.J[ )- |^73| ), which 
belongs to k2.2) . and, in addition, 

Each weak solution satisfies 



(4.5) (5tu(0,w)(^i(f^)).^^i(j^)+ J \'A,{u{t))-\'w-u{t){VIC*u{t))-\'wdx = 0, 

for all w e {Q), for almost all t G [0, T] , where (y) := A (y) + ey, £ > 0. 
Remark 4.2. The weak solution of (|4.ip - (14.21) satisfies the energy identity 



£eiu{t))+ / u{s) V$^(m(s)) - V/C*u(s) dxds^S^iO), 
Jo Jn 

for all t >Q, see Remark l2.6l The critical mass Mc > in the statement of Theorem 
14.11 is estimated in [171 Proposition 2 and 3]. 

In view of estimate (|3.36p . the following result can be proven for (|4.ip - (|4.3p . 



Lemma 4.3. Let the assumptions of Theorem \4-4\ satisfied. Then, for every 
T > 0, there exists a constant Cm,s.t > such that 

(4-6) sup ||u(t)||^i(j^) + ||atu||^2([t_t+i]xn) <CM.e,r- 

Moreover, for any bounded set B C Zun , there exists a time = {B) > such 
that S, (t) B C {n) , for all t>t^. 

Proof It suffices to show (j^B]) with t = t+ + 1 (see ([04| - (|05| ). where t+ > 
is the time given in Theorem 13.51 (which still applies for A^ (y) ~ A(y) + ey) . Set 
A^ {y) — A^ (s) ds. Testing equation ()4.5p with w — A^{u (t)), and integrating 
over SI, we deduce 

d ,^ /■ ,_ 2 



(4.7) - {A, (u (t)) , 1)^,(^) + J^^ \VA, (u m' dx 

u (t) (V/C * u {t)) ■ \/A, {u (t)) dx 



n 

<\j^ |VA, {u {t)^ dx + C \\u IIV/C * u 

Integrating this inequality from t to t + 1, and using the fact that u is bounded 
according to p.36p , we obtain 

(4.8) [ \VA{u{s))\^dx<C, 

Jf Jn 

for all t > t+ (M) , for some positive constant C independent of time and the initial 
data. 

In order to rigorously prove ()4.6p . one must proceed as for the problem (|2.6p . 
More precisely, recalling that problem (|4.ip - ()4.3p is uniformly parabolic, one has to 
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employ another regularization scheme in which is approximated by a sequence 
of functions (Ae)^ e C°° (R+) , the data u\t=Q = uoe G (fl) n (IT) is such 
that uoe — > uq in the L°°-metric, and /C is replaced by a sequence of smooth kernels 
/Ce G C°° such that /Ce — >■ /C in W^^^ (W^^ , as e 0. The procedure ensures that the 
approximate solutions (u^)^ are smooth enough so that all the computations below 
can be performed rigorously. Thus, in what follows we shall again proceed formally 
(it will be easy to see that all the constants in the estimates are independent of 
e — > 0+). To this end, testing equation (|4.5|) with w — dt {A^u {t)) , and using the 
fact that A^ {y) > £ > 0, we deduce 

(4.9) 2e [ \dtu\^dx+^ [ \VA,{u)fdx:^2 [ u{WK * u) ■ WdtA, {u) dx. 
Since 

/ u (V/C * u) ■ VdtA, {u)dx^ ^ I u (V/C * u) ■ VA^ [u) dx 
Ja dt Jq 

dtu (V/C * u) ■ VAg (u) dx 



- u (V/C * dtu) ■ VAe (u) dx, 
Jn 



relation (14.91) together with basic Holder and Young inequalities imply 

(4.10) e [ \dtufdx + 4- f \VA,{u)fdx 

Jn dt Jq 

<4 / 2u{VIC^u) ■VA,{u)dx + c [ \VA,{u)fdx, 
dt Jn Jn 

for some positive constant c > independent of time and the initial data. Next, 

we multiply (|4J0l) by e"^''^^^ for s e {t,t + l), to obtain 

(4.11) Ae^(*--) / \VA,{u)fdx<e<'-'^-^ f 2u (V/C * u) • VA, (m) dx. 
ds Jn ds Jn 

Integrating (|4.11l) between s and t + 1 gives 

(4.12) /" \VA,{u{t + l))fdx~e''^'-''^ [ \V A, {u is))f dx 
Jn Jn 

< / e'=(t-^)_:i / 2u{T){WIC*u{T))-WA^{u{T))dxdT. 
Js dr Jq 

Notice that we can split the integral on the right-hand side of (|4.12p as follows: 

/ '2u{T){VK*u{T))-VA,{u{T))dxdT 



dr 

e=(*-^) / 2u(T)(V/C*w(T))-VAe(u(T))dx|*+i 



t+1 

(-c)e'=(*"^) / 2u{t){VIC*u{t)) ■VA,{u{T))dxdT 



=■■ h + It 



Next, 



'^t*-^) / 2u{t){VIC*u{t)) ■\'A,{u{T))dx\ 



t+1 
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= e-'^ / 2u{t+ 1) (V/C * w (t + 1)) • VAe (u (< + 1)) dx 
_gc{t-s) f 2u{s){\/JC*u{s)) ■\/A,{u{s))dx, 
which can be further bounded, exploiting standard Holder and Young inequalities, 

by 

(4.13) ^6 - V X ^ ^ 

Moreover, we have 



rt+l 

(4.14) /7 = / ce'^f*-^) / 2u(T)(VA:*w(r)) • VA, (u(T))da;dr 



t+i 



< C / ( / |u(r)|^da;+ / \V A, {u {t))^ dx) dr. 

Thus, on account of (|4.13l) - (|4.14l) and estimate (|4.19p . inequality (|4.12p becomes 

(4.15) e"" / \VA,iu{t+l))fdx~e''^*-'^ [ \V A, {u {s))f dx 

Jn Jn 

< ^ / \VA,iu{t + l))fdx+ f \VA,iuis))fdx 
^ Jn Jn 

+ C ( [ \u{Tfdx+ [ \WA,{u{T))\^dx]dT + C, 



n Jn 



for all t> t+. Therefore, 



^— [ \VA,iu{t+l))fdx 
^ Jn 



<e''(.t-s) / \\jA,{u{s))\^ dx+ / \VA,{u{s))\^ dx 
Jn Jn 

+ C ( [ \u{T)\^dx+ [ \WA,{u{T))\^dx]dT + C. 
Js \Jn Jn J 

Integrating this inequality from t to t + 1 with respect to s, and recalling (|4.19p we 
obtain 



(4.16) 

^ / \VA,iuit + l))fdx 
^ Jn 



n 

< / e<'-''^ / \VA,{uis))f dxds+ / / \VA,iu{s))f dx 
Jt Jn Jt Jn 



t+l rt+l 



C I I [I \u{Tfdx+ / \\/A,{u{T))\^dx\ dr + C 



It Js \Jn Jn 
l-t+i !■ ^t+i 

<C \\7Ae{u{s))fdx + C (r-t) \\/Ae{u{T))fdxdT + C. 
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By virtue of gH), (|4J6l) yields 

(4.17) / |VAe {u {t + l))\^dx <C, yt> t+. 

Since [y] > e, for all y, and integrating (I4.10[) over (t, t + 1) once more, (|4.17p 
entails the desired estimate (14.61) . The proof is finished. □ 



As in Section 3, we can prove the following result for (I4.ip - (|4.3p . 



Theorem 4.4. Let the assumptions of Theorem \4-l\ be satisfied. The dynamical 
system [S^ (t) ,Zjj^) ({S,; (i)}j>o defined as in {S.^b^ ) possesses a global attractor 
A — As,M in the sense of Definition \3.6\ such that A is globally bounded in C" {p,)r\ 
(n) , for some a e (0, 1) . 

Proof. As in the proof of Theorem 13.71 we must check that (t) is closed and 
that it admits a compact absorbing set in Z^n- 

Step 1 (Closedness of Se). The proof of this step is essentially the same as in 
Lemma 12.71 where everywhere in the estimates we must replace A by the function 
Agr. This only affects the estimate for (|2.20p . which now reads 



(4.18) h = - (Ae (wi) - A^ (m2)) (wi - U2) dx + M12 / A, (ui) - A, (ua) dx 



n 



-e \\ui - W2|li2(o) - / - ^ ("2)) {ui - U2) dx 



+ £ {M12) + *^i2 / A {ui) - A (lis) dx 
Jn 

< -e \\m - U2\\l.^n) +C{l + e) {A'h - I^hf . 

Thus, the same inequality in (|2.15|) is valid for any two weak solutions ui (t) , U2 (t) 
of problem (|4.1|) - (j4.3|) corresponding to the initial data uio,W20- 

Step 2 (Smoothing effect). This step requires only minor modifications in the 
proof of Theorem 13.51 Indeed, Lemma 13.31 also applies to the function A^ and 
the results in, e.g., [20] Corollary 4.2] can be still applied to obtain global Holder 
continuity of the weak solutions. In particular, each weak solution of (|4.ip - (14.3p 
satisfies the L^-C^DH^ smoothing property, and there exists a time t+ — t+ {M) > 
such that 

(4-19) sup ||M|lc»/2,c([t,t+i]xn)nL~([t+;oo);_f/i(n)) ^ C'm, 



+ 



for some positive constant Cm, independent of time and the initial data. Hence, 
a compact in Zdn absorbing set like in p.50p can be sought, which is enough to 
apply [461 Corollary 6] once again. The proof is finished. □ 

Next, taking advantage of the asymptotic smoothness of (t) we can also show 
that (14.ip - (j4.3p has a gradient structure. To this end, we define the w-limit set of 
a trajectory u (t) of (|4.11) - (|4.3p . starting from 

uo e Zdn = {uo e L°° (fl) : uo > 0, (uq) < M}, 

as follows: 

(4.20) io {uq) = |u* e Zdi^ : 3i„ — > od such that lim \\u (f„) — u*|lioo(Q) = o| . 
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We prove that ui (uq) consists of stationary solutions, satisfying the system 

{div ^u»V (u*) — /C * = 0, in fi, 

(li,) - -'C * "*) ■ ^ 0, onr, 
jf^T "* (^) '^^ ("o> , 
where we recall that $g (y) = (y) /y such that $^ (1) = $e (0) — (see Remark 

The following proposition justifies to call 
(4.22) £,{u{t))= [ ^,{u{t))dx~~^ [ [ u{x,t)lC{x~y)u{y,t)dxdy 



an energy functional for (|4.ip - ()4.3p . 



n Jo 



Proposition 4.5. Let u(t) = S^{t)uo, uq £ Z^n, be the (unique) global so- 
lution of the non- degenerate aggregation equation |^._?[ )- |^3[ ). Then the following 
assertions are true: 

(i) The function (u (t)) is differentiate a.e. on [r, oo), for every r > 0, and 

d f 

(4.23) —£,{u{t))=^- u{t)\7q>'^(u{t))-\/IC*u{t) dx, 
dt 

for a.e. t > 0. 

(a) The function £^ [u (t)) is nonincreasing , and there exists a positive constant 
C,, depending only on JC, Q, e and \\u\\j^oo(^^-^ such that 

(4.24) Se{u{t))>-C,, 

for a.e. t > 0. Moreover, limt_i.oo £e (u (t)) = infoo £s {u (t)) = ^e.oa G R exists. 
(Hi) For any uq £ Zun , w [uq) is a compact, connected invariant set, and every 
G w (uq) is a solution of the stationary problem \^.21\ l. 

(iv) Every £ uj (uq) is a critical point of £^ [u (t)) in ^.22\ l, i.e, £^ (u*) = 0. 

Proof. First, note that by the assumption (H2) and (|4.4p . we have ^ {y) ^ A (y) ^ 
y™ with TO > TO*, and by definition, $e (y) — $ {y)+s {y In (y) — y) , for every e > 0. 
The first part of assertion (i) follows from the fact that $e G and the regularity 
of the bounded solution u (t) on the intervals [r,oo), for every r > (indeed, u G 
C"/^ (r, oo; C" (fl))). The second part is a consequence of Remark li?^ By 
£e {u (t)) is nonincreasing on Z^n. The elementary inequality y In (y) — y> —1, for 
all y > 0, and the assumption (H2) on A, yields on Zdn that <I>e {u (t)) > — Cg, for 
some constant > depending only on the physical parameters of the problem 
and the L°°-bound of u. Rewriting the energy £g (t) in the following equivalent 
form 



1 



£,{u{t))= I <^>,{u{t))dx + ^ I I lC{x-y){u{x,t)~u{y,t)f dxdy 

-7; I a{x) {u{x,t)f dx, 
^ Jn 

where 

a{x)^ [ JC{x~y)dyeC (Ti) , 

the second part of assertion (ii) is also immediate. First, by Theorem 14.11 problem 
()4.ip - (|4.3p defines a nonlinear Co-semigroup {t) : Zdn Zjjn, u {t) = {t) uq 
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with uq E Zd n ■ Second, by the resuhs of Theorems 13.51 and 14.31 we know that for 
any uo G Zdn, there is i+ > such that [Jt>t^Se {t)uo is bounded in {il) n 
C" (f2) , and hence relatively compact in Zdn (when endowed with the metric 
topology of L°° {fl}). Third, it can be seen from (i) that the function {u{t)) : 
J'DN — > M is a Lyapunov function on Tdn C Zjjn, for any (positively invariant) 
subset J'DN C (il) nC" (ft). In particular, by Theorem l4.4l we can take J-dn = 
-^e,M, where Ae,M is the global attractor for {S^ (t) , Zdn) ■ Moreover, satisfies: 
if for t > 0, £e(S'e(t)u*) — £e (u,) , then it, is an equilibrium point of Se{t). 
In conclusion, by [IHl Chapter 10, Definition 10.1], {S^ (t) ,Ae,M) is a gradient 
system. Thus, by HH Propositions 10.3 and 10.12] we immediately conclude (iii), 
i.e., io (uq) is a nonempty, compact, connected invariant set, and lu (uq) consists 
only of stationary solutions. The final part of (ii), inft>o Se {u (t)) — (u^) — £e,ca 
is satisfied owing to and (iii). 

Finally, for (iv) we observe that if e a; (uq) is a solution of (I4.2ip . then for 
any i/> e 77^ {Q) n L°° [Q) it follows from ^^jTi\ that 

= y div ^u, V (u,) — /C * u,^ ^ ipdx — J V (^^^ (u^) — K. * m*^ • itipda 

(VAe (u,) — u,V/C * u,) • Vxpdx, 
which, by straightforward computations, is just the following: 

(w, +<5V')|5=o = 0' 

i.e., is also a critical point of £^ in Zdn- In fact, we easily see that the statements 
(iii) and (iv) are equivalent to each other. The proof is complete. □ 

Next, we show that uj (uq) has a positive bound from below depending only on 
the physical parameters of the problem and uq. 

Proposition 4.6. Let u (t) — (t) uq, uq G Zdn, be the unique solution of J[ )- 
li4-3^ such that (uq) > 0. Suppose that E uj (uq)- Then, there exists a constant 
u > 0, depending only on uq, Q, K. and £ > 0, such that (x) > u > 0, for all 
X € il. 

Proof. The proof follows a similar argument used in '26' Section 2] , [54j Proposition 
3.3]. Indeed, owing to the homogeneous Neumann boundary condition of (|4.2ip . 
the first equation of (|4.2ip also reads 

J u*V (^<p'^ (u*) - /C * u,^ • \/^pdx = 0, 

for any ip E (fi), whence V($g (w*) — /C * u,) = in each connected component 
of the open set where x E such that {x) > 0. Therefore, (u,) — /C * is 
constant on each such connected component. If (uq) > 0, by the last equation of 
(|4.2ip . there exists G f2 such that (xq) > 0. The previous statement implies 
that there exists 7 e M such that 

(4.25) 7 = $e (u, (x)) ~ (/C * u*) (x) = (u* (x)) + £ In (u, (x)) - (/C * u,) (x) , 

for X e il (xo) , which is the same connected component of {x g f2 : it, (x) > 0} as 
Xq. Since u, G w (mq) is bounded and /C G Wj^'J^ (K'*), we observe from (|4.25p that 
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u* (x) satisfies the inequality 

i('7-||^*u.lL = (n)-||*'KW)|l _ 

u* (x) > e V ^ ' II iii°=(n)y > 0, 

for all X € (xq) . In particular, this yields that fi (xq) = f2 and the claim is 
proved. □ 

Remark 4.7. Some results on properties of the steady states for the aggregation 
equation with nonlinear diffusion in the case f2 = R'' can be found ^ & |9j . 

As a result of the proof of Proposition 14.51 we can now conclude the following 

Theorem 4.8. Let the assumptions of Theorem \4-4\ be satisfied. The global at- 
tractor A = Ae,M of problem consists entirely of unstable manifolds of 

the equilibria, which are bounded solutions of \4-.21^ . 

Under an additional assumption on the kernel which ensures that K, is reasonably 
smooth at the origin, we can show that Ae.M is also finite dimensional. 

Theorem 4.9. Let the assumptions of Theorem \4-.l\ be satisfied. Assume 

(4.26) {D^IC)lB,^o)eL^{R'i), d>2. 

The global attractor Ae.M of ^J7J^-^J^ has finite fractal dimension: 

dimp {Ae,Ai,H^^ (n) n LP (n)) < CeM < oo, 

for any 1 < p < od. 

Remark 4.10. Recall that fractal dimension of a compact set y, dimp {y,X) is 
defined as 

dimF(3^,^)= hm — — — — , 

5^0+ log2 (1/(5) 

where Ng {y,X) is the minimal number of balls Bg that can be used to cover the 
compact set y in the metric of X. 



The statement of Theorem l4.9l is in fact a consequence of a much stronger result 
which states that (|4.ip - (l4.3p admits an exponential attractor Aie^M provided that 
(|4.26p is also satisfied. The precise statement is given by the following. 

Theorem 4.11. Let the assumptions of Theorem \4-.9\ be satisfied. Assume \4.26^ . 
For every fixed e > 0, there exists an exponential attractor M = Ai^.M bounded in 
C" (ri) for the dynamical system {S^ {t) ^Zdm) which satisfies the following prop- 
erties: 

(i) Semi-invariance: (t) C M, for every t>0. 

(ii) Exponential attraction: For every bounded subset B C -Zdtv, 

(4.27) distHi-f^n)nL^(^n){Se[t)BM)<Ce-'^\ > 0, 

for some positive constants C — C (e, M) and k, for any 1 < p < oo. 
(Hi) Finite dimensionality: 

(4.28) dimj. {M,H^- [VL] n L^ {n)) < C^m < oo, 

for any 1 < p < oo. The constants C,Cs;,m,I'^ can be computed explicitly in terms 
of the physical parameters of the problem. 



MODEL FOR BIOLOGICAL AGGREGATIONS 



27 



Here, and everywhere else, we denote H^^ {il) H^^^ {ft), for any S £ (0, 1]. 
Since the global attractor Ae^ai is always contained in Me.M, the above theo- 
rem immediately implies that the fractal dimension of the global attractor Ae^m is 
also finite. It is worth to recall that, in the global attractors theory, it is usually 
extremely difficult (if not impossible) to estimate and to express the rate of conver- 
gence in p.47|) in terms of the physical parameters of the system considered. This 
constitutes the main drawback of the theory. Simple examples show that the rate 
of convergence in p.47p can be arbitrarily slow and non-uniform with respect to the 
parameters of the system considered. As a consequence, the global attractor be- 
comes sensitive to small perturbations and, in some sense, cannot even be observed 
in experiments. The concept of exponential attractor overcomes this difficulty (see, 
e.g., the survey article [4F). Indeed, in contrast to the global attractors theory, the 
constants C, Ce,M, k in (|4.27p - ()4.28p can be explicitly found in terms of the physical 
parameters. 

We report for the reader's convenience the following abstract result on the 
existence of exponential attractors [23 Proposition 4.1] (cf. also [28l Proposition 
2.17]) which will be used in the proof of Theorem 14. Ill 

Proposition 4.12. Let "HjV.Vi be Banach spaces such that the embedding Vi C V 
is compact. Let B be a closed bounded subset of % and let S : B ^ B be a map. 
Assume also that there exists a uniformly Lipschitz continuous map T : _B — > Vi, 
i.e., 

(4.29) IITOi-Tballv, <L||5i-fe2|l„, Vbi^b^ E B, 
for some L > 0, such that 

(4.30) \\Sbi-Sb2\\n<e\\bi-b2\\.H+K\\Tbi-Tb2\\v, ^biM^B, 

for some 6 < ^ and K > 0. Then, there exists a (discrete) exponential attractor 
Add (1 B of the semigroup {S{n) §",n G Z+} with discrete time in the phase 
space %. 

We delay the proof of Theorem 14.111 until the end of the section. The idea is 
to verify (|4.29p - (|4.30p for a suitable choice of maps. We begin by showing that the 
semigroup {t) is strongly (Lipschitz) continuous with respect to the (iJ-^ (17)) - 
metric. 

Proposition 4.13. Let Ui, i = 1,2, be a pair of weak solutions according to the 
assumptions of Theorem \4.9\ Then the following estimate holds: 

(4.31) {t) ~ U2 (i)ll(Hi(j^))- + e / hi (s) - "2 (5)11^2(0) ds 

•J 

< hi (0)-^.2 (0)||^^,(^)).e«*, 

for all t>Q, for some positive constants n, C which depend on e > and K, but are 
independent of Ui (0) . 

Proof. Following Lemma 12.71 we have that u := ui — U2 and ri[t) :~ \W{t)\\^f^fjiY 
satisfies the problem 

\j^n (t) = - {dtu (t) , (t)) =ii + i2 + 13, 
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with Iijh, h given in ()2.19p . The integral Ii is estimated in ()4.18p , whereas for 
12,1?, we have (12.2311 . The estimate for I2 can be improved in ()2.23p . using (|4.26p . 
Young's inequahty for convolutions and Lemma |2. 11 (b), as follows: 



I2<C j \D^JC * ui I |V0|^ dx<C \\D^JC * ui ll^^(j^) II V0II 



2 



|2 



< C [\\D ^||ii(Bi(o)) II"IIl~(o) + 11^ ^ILp(r^\Bi(o)) ll"llLp/(p-i)(n) J II ^ vwL^n) 
Thus, we get 

d 
di 

which yields the desired inequality (|4.31l) by application of the Gronwall's inequality. 

□ 



(4.32) f\\uml^,y+2e\\u{t)\\l.,<C\\u{t)\f^^, 



Remark 4.14. A crucial point in the proof of Theorem 14.111 is that we need the 
global Lipschitz continuity of 5^ (t) in the norm of (if^ (17)) . The assumption 
(j4.26p plays an essential role with respect to this issue (see (j4.32p ). While Newtonian 
potentials do not satisfy (|4.26p . in population dynamics the non-local effects are 
generally modelled with smooth, fast-decaying kernels K, which obey (|4.26l) . see 
e.g., Pl. 

The step needed to establish the existence of an exponential attractor is the 
validity of so-called smoothing property for the difference of two solutions of (14. ip - 
(|4.3p . In the present case, such a property is a consequence of the following two 
lemmas. The first result establishes that the semigroup S'^ (t) is some kind of 
contraction map, up to the term ||mi — '"2|Il2([o t]-{H^)*)- 

Lemma 4.15. Let the assumptions of Proposition \4-l^ hold. Then, for every 

t > 0, we have: 

(4.33) 



||ui (t) - 1*2 (t)ll(^i). < e " ||mi (0) - U2 (0)11(^1). +C£,Af J ||ui (s) - U2 (s)ll(^i). ds, 

for some positive constants C^.m, k which depend on e > 0, Q and K,. 

Proof. Recall that u := ui — U2. Combining (j4.32p together with Poincare's in- 
equality 

\a-"^ {u~{u))\^^^ + {uf<Cn\\u\\l., 
we deduce from (|4.32l) the following inequality: 

^ll"WII(Hi)* +|^h(i)ll(//i)- <C||u(t)||^^i)., 
for all t > 0. Thus, Gronwall's inequality entails the desired estimate (|4.33p . □ 

We now need some compactness for the term ||ui — U2|li2(-[Q on the 

right-hand side of (|4.33p . This is given by 

Lemma 4.16. Let the assumptions of Proposition hold. Then, for every 

t >Q, the following estimate holds: 

(4.34) ll^t^i -^t"2|li2([o_t].i)(A„)') ||ui (s) - U2 (s)||i2(f^) ds 
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<C,Me'''\\ui (0)-U2 (0)11(^1). , 
where Ce,M > and k > also depend on e, Q and IC. 

Proof. The second term on the left-hand side of (|4.34p can be easily controlled by 
(|4.3ip . Thus we only need to estimate the time derivative. Recall that each dtUi, 
i = 1,2, satisfies (|4.5p . Furthermore, in light of Theorem 14. 11 recall that we have 

(4.35) sup||u,(i)||ioo(f2) <Ce,A/, * = 1,2. 

t>0 

Thus, for any test function w G D{An), using the weak formulation (j4.5p . for 
dtu :— dtUi — dtU2, there holds 

{dtu {t),w) = h+h, 

where 

h ■.= -{V{A,{ui)^A, iu2)),Vw), 

h ■■= {U (V/C * Ml) - U2 (V/C * u) , Vw) . 

First, for every w Cz D (Aj^) we have 

h = {A^ (ui) - A^ (u2) , Aw) < Q ||ui - U2|li2 ||Aw||^2 

< C||wi -U2||^2 \\w\\d(Aj^) ■ 

On the other hand, it is easy to show, on account of (|4.35p . that 

h <C \\ui - U2||i2 ||Vu)||^2 . 

These estimates together with (|4.3ip gives the desired estimate on the time deriv- 
ative in (gJll). □ 

We now show that the semigroup Sg (t) is actually uniformly Holder continuous 
in the H^^ L^-norm with respect to the initial data. 

Lemma 4.17. Let Ui (t) = S {t) Ui (0), with Ui(0) G Zdn- Then, for any 1 < p < 
oo, the following estimate is valid: 

(4.36) \\u, (t) - U2 (i)||^i-(o)nLp(0) < Ce,Me'^* \\u, (0) - mijn^r ' 

for all t > t^,, where the constants Ce,j\/ > 0, k > and J — J (p) < 1 are 
independent of the initial data and time. 

Proof. Using the interpolation [H^, (i?^)*]i/2,2 — L'^, we deduce from estimates 
(|43T|1 and (jM]) that 

(4.37) hi {t) - U2 (t)|U2(j,) < CeMe'^' hi (0) - m\\'my , 

for all t > = t+ + l. By interpolation in the spaces L°° C C L^, for 2 < p < cx), 
fji ^ jji-s ^ (iJi)*, S G (0, 1), the estimate (|436ll also holds for the difference of 
solutions u — ui — U2. □ 

The last ingredient we need is the uniform Holder continuity of f i— )■ Sf;(t)uo in 
the n L^-norm, namely. 
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Lemma 4.18. Let the assumptions of Proposition \4-13\ be satisfied. Consider 
u (t) = Se {t) uq with uq e Zdn ■ The following estimate holds 

(4.38) \\uit)^u{s)\\^,^^^^^^,^^^<C,,M\t~s\\ yt,s>t,, 

where 7 = 7 (p) < 1 and the positive constant Ce,M is independent of initial data, 
u and t, s. 

Proof According to and (gH), and recalling that K. G W^^^ {R'^) , the fol- 

lowing bound holds: 

sup \\AA,u (i) - V • (u (V/C * u (t)))|| f^i,. < Ce,M- 
t>t. 

Consequently, by comparison in (14. 5p . we have that 

sup \\dtU {t)\\,^iY < Ce,M, 

which entails 

(4.39) \\u{t)-u{s)\\^J^^y <CeM\t-s\, Vi,s>t,. 

Estimate (|4.38p now follows from (|4.39p . and standard interpolation inequalities, 
as in the proof of Lemma 14.171 □ 

Proof of Theorem \4-ll\ In order to apply Proposition 14.121 it is sufficient 
to verify the existence of an exponential attractor for the restriction of S{t) on 
some properly chosen semi- invariant absorbing set in Zjjn. Recall that, by ()4.19p 
and Lemma 14.31 the ball Bq := B^^^-^-^^^i^^s^ {Cs^m) will be absorbing for (t), 
provided that C^^m > is sufficiently large. Since we want this ball to be semi- 
invariant with respect to the semigroup, we push it forward by the semigroup, by 
defining first the set Bi = [Ut>oSe {t) i?o](//i)*, where denotes closure in the 

space (i?^ i^))* ■ Then set M = S (1) Bi. Thus, B is a semi-invariant compact (for 
the metric of [H^ i^))*) subset of the phase space Zdn- On the other hand, due 
to the results proven in this section, we have 

(4-40) sup(||ii(t)||p„(^)^^i(,j) + \\A,u{t)\\^,^,,^ + 119*7/(011(^1)-) < Ce.M, 

for every trajectory u originating from uq G B, for some positive constant C^^m 
which is independent of the choice of uo S B. We can now apply the abstract result 
above to the map S = S'^ (T) and "H = (if^ i^))* , for a fixed T > such that 
e''^'^ < 5, where k > is the same as in Lemma 14.151 To this end, we introduce 
the functional spaces 
(4.41) 

Vi ([0, T] ; (n)) n ([0, T] ; DiA^y) , V ([0, T] ; {H^ (17))*) , 

and note that Vi is compactly embedded into V. Finally, we introduce the operator 
T : B ^ Vi, by Tuo := u e Vi, where u solves with u (0) = uq e M. 

We claim that the maps S, T, the spaces H,V,Vi thus defined satisfy all the as- 
sumptions of Proposition 14.121 Indeed, the global Lipschitz continuity (|4.291) of T 
is an immediate corollary of Lemma [4. 161 and estimate (|4.30p follows from estimate 
(|4.33p . Therefore, due to Proposition 14.121 the semigroup S{n) = Se (nT) gener- 
ated by the iterations of the operator § : B — >■ B possesses a (discrete) exponential 
attractor Md in B endowed by the topology of [H^ (fl)) . In order to construct 
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the exponential attractor £ for the semigroup Se (t) with continuous time, we note 
that, due to Lemma 14.131 this semigroup is Lipschitz continuous with respect to 
the initial data in the topology of (H^ i^))* ■ Moreover, by (|4.36p and (|4.38p the 
map (t, uq) i-^ Ss (t) uq is also uniformly Holder continuous on [0, T] x B, where B 
is endowed with the metric topology of (i?^ (^^))*- Hence, the desired exponential 
attractor A4 for the continuous semigroup S^it) can be obtained by the standard 
formula 

(4.42) \J SAt)Md. 

*e[o,T] 

In order to finish the proof of the theorem, we only need to verify that M defined as 
above will be the exponential attractor for (t) restricted to B not only with respect 
to the (i?^ (r2))*-metric, but also in with respect to a stronger metric. This is an 
immediate corollary of the following facts: B is bounded in C" (fl) n (f2) and 
standard interpolation inequalities between the following spaces: L°° C LP C L^, 
iJi C H^-^ C {H^Y , for 2 < p < oo and ^ e (0, 1]. Theorem iH] is now proved. 

5. Convergence to steady states 

In this section, we show that any global-in-time bounded solution to the model 
(|4.ip - (l4.3p converges to a single equilibrium of (|4.2ip as time tends to infinity. The 
proof of the main result is based on a suitable version of the Lojasiewicz-Simon 
theorem and Propositions 14.51 14.61 The question of such convergence is usually a 
delicate matter since it is well known that the topology of the set of stationary 
solutions of (|4.2ip can be non-trivial. In particular, there may be a continuum of 
stationary solutions for (I4.2ip even in the simplest cases, for instance when is a 
disk, IC is either a Newtonian or Bessel potential and A (y) = y (see [39l [50]). 

The main result of this section reads as follows. 

Theorem 5.1. Let the assumptions of Theorem \4-l\ be satisfied. Assume that $ is 
a real analytic function on R-|_. For any uq G Zdj^ with (uq) > 0, the corresponding 
positive solution u (t) = (t) uq to the non-degenerate aggregation equations (^T7p" 
Iji4-3{ l converges to a single stationary state of 14-21^ in the sense that 

(5.1) ^lim ||u(t)-u,||^p(j^j =0, 

for any p > 1. Moreover, there exist constants C > 0, p — p{p,u,,) G (0, 1) such 
that the convergence rate holds: 

(5.2) ||7.(t)-u.||^,(,,) <C(l + t)-^ 
for all t > 0. 

Owing to ()4.20p . Proposition 14.51 (iii) . Proposition 14.61 ^nd the regularity prop- 
erties of M e C"/2 ((0;oo);C" (n)), we may then assume without loss of generality 
that 

(5.3) inf u{t,x)>u> 0, for ah t > 0. 

We employ a generalized version of the Lojasiewicz-Simon theorem proved in 
[5T1 Theorem 6] (cf. also [IHl Lemma 2.20]). The version that applies to our case 
is formulated in the following. 
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Lemma 5.2. Under the assumptions of Theorem \ 5.1\ there exist constants 6 € 
(0, C > 0, (5 > such that the following inequality holds: 



(5.4) \£s (u) - £e < C ll/i - ifAWmn) > 

for all u G L°° {fl)r]H^ (O) provided that ||u — |1^2(j7) < S. Here fJ, — (u) denotes 
$g (u) — /C * u. 

Proof. We will now apply the abstract result [311 Theorem 6] to the energy func- 
tional £g (u), which according to (j4.22l) is the sum of entropy and an interface energy 
term. In contrast to this feature, we shall split £5 (ip) into the sum of a convex (en- 
tropy) functional S : {fl) — > M U {00}, with a suitable effective domain, and a 
non-local interaction functional 5* : (fi) — M. Let u > be the lower bound 
from (|5.3I) and recall that e > 0. We define the lower-semicontinuous and strongly 
convex functional E = by 

^ ($e (") - e In (u) u) dx, iiueL°°{n), u> 0, 
-|-oo, otherwise, 

with closed effective domain dom(I]) = ZdnDH^ (fl), and the quadratic functional 
* = *e : ^ given by 

* (u) / [u{JC*u) - 2£ln(u)u] dx. 

2 Jo 

We have that E is Frechet differentiable on any open subset U of 

Um W € L°-{n) : K^-)! < M; u<^{x)< Cm} , 

(where Cm > is such that || V'IIloc(o) < Cm, since 1/; is bounded) with Frechet 
derivative DE : U — ?> L°° (fi) having the form 



E(u) 



(ra(7.),0= / (<(«)-eln(2i)) -edx, 



for all u £ U and ^ G (f2). The analyticity of I?E as a mapping on (f2) is 
standard owing to the analyticity of $ (see, e.g., [311 Remark 3]). Moreover, due 
to assumptions on A in Theorem 14.11 '^e have $ (y) ^ A (y) ^ y™^^ and recalling 
that $g (y) = $' (y) -I- £ In (y), one has $' (u) > and 

(DE (ui) - DE (M2) , wi - M2) > * (u) 11^1 - W2||^2(Q) , 

for all Mi,M2 G U, and 

IIDE (ui) - DE (w2)||(i2(o))- < 7 hi - W2|Il2(j^) , 

for some positive constant 7 = 7(Cm,w, £)■ Moreover, computing the second 
Frechet derivative D^I] of E, 



(D'^iu)^,,^,)^ [ ^':{u)^,-^^dx 



yields that D^T, G C {L°° (fi) , L°° (il)) is an isomorphism for every u £ U, owing to 
the fact that (u) = (u) /w > ^ = ^ {Cm,u} > 0. Concerning the (quadratic) 
function ^P, we see that 

*(") = ^ (-/C *u,u)^2(s:2) + (eln(u),w)^2(f^) , Vu G (fi) . 
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We recall that the linear operator i/j i— > /C * is self-adjoint and compact from 
i2 to itself and is also compact from L°°{n) to C°{Ti) (since K. € W^;^). On 
the other hand, we also have the following (orthogonal) sum decomposition of 
(fl) = Ll {n)®Hi, where 

Ll {VL) -.^ {u e {n) : (u) =0} , Hi := {u G (Q) : u = const.} 

Then, the annihilator of Lq (fJ) is the one-dimensional subspace 

Ll \ch e {L^ (fl))* : c e ]r| 

of constant functions, where h g (L^ (^))* — (^) is given by {h, u) = udx, 
u £ L'^ (fi) . Hence, the hypotheses of |3ll Theorem 6] are satisfied and the sum 

= S + * : ^ K u {oo} 

is a well defined, bounded from below functional with nonempty, closed, and convex 
effective domain dom{£^) =dom(S) . Unravelling notation in |3H Theorem 6], and 
observing that the Frechet derivative 

(u) — {u) ~ JC * u —: /i, 

we have 

(u) -£,{u.)\'-' <C inf ( (u) - fij^^n) ■■ ^ L^} 

IL ^ ( i * ) V J 

from which (15.41) follows. □ 



We are now in a position to complete the proof of Theorem 15.11 

Proof of Theorem \5.1\ We argue as in the proof of [2H1 Theorem 2.21]. First, we 
note that by virtue of the regularity results proven in the previous section (see, 
e.g., (|Tig)) ). ah e uj[uo\ are bounded in C" (H) n (f2). Besides, recaUing 
Proposition 14. 51 - (ii) . we have 

£e {u (t)) £e,oo, as i OO, 

and the limit energy Ee,oo is the same for every steady-state solution m* £ a; [uo]. 
Moreover, we can integrate ()4.23p over (i, oo) to get 



(5.5) u / V^,_{u{t))-V]C*u{t) 



dxds 

< / u{t) {u (t)) - V/C * w (t) dxds 
Jt Jn 

= £e {u (t)) - £e,oo = £e {u (t)) ~ £e (u*) . 

By virtue of Lemma [5.21 and recalling that ^ [t) ~ $^ [u [t)) — V/C *u{t), we have 

(5.6) % {u {t)) - {u,)\^-' <C\\^i{t)- {fi (t))|L.(,,) < C IIVa. mmn) 
exploiting Poincare's inequality, provided that 

(5.7) ||m-u,||^2(s:^) <(5. 
This, combined with the previous identity, yields 



/oo 
||VM(s)||i.(,,) dS<C\\Wpi{t)\\]^y 
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for all < > 0, for as long as (|5.7p holds. Note that, in general, the quantities 9, C 
and 5 above may depend on > and e > 0. Finally, let us set 

= U {X : Z is an open interval on which (|5.7p holds} . 

Clearly, W is nonempty since G t^[uo]. We can now use (|5.8p . the fact that 
|lV/i(t)||^2(Q) G L2(0,oo), and exploit |M1 Lemma 5.1] (with a = 2(1-6*)) to 
deduce that || (•)||^2(q) G L'^{W) and 

(5.9) / ||V/i (s)||^2(f^x ds < C (u*, u) < cx). 

Consequently, using the bound (|5.9p and the main equation (|4.ip , which also reads 
dtu (t) — div {u (t) V II (t)), wc obtain 



(5.10) j^\\dtu{s)\\^^,^^^yds<C 

In order to finish the proof of the convergence result in (15.11) it suffices to show 
that it holds in L^-norm. Indeed, in this case (15.11) will become an immediate 
consequence of the L^-(C" n i?^) smoothing property of the bounded solutions 
u it) and all G w [uq] • We claim that we can find a sufficiently large time r > 
such that (t, oo) C W . To this end, recalling (|5.5p and the above bounds, we also 

have that dt^ G L^{Q,oo] {H^ (^2))*), V^i G ^^(o^oo; (L^ {^)Y) and, furthermore, 

for any fc > there exists a time = {k) > such that 

(5.11) 

ll^t^l!Li(W'n(t,,oo);(J?i)*) - ^' ll^t"llL2((t,^oo);(Hi)*) - ^' 1 1 II ((t, _oo) ;(L2)d) < k. 

Next, observe that by the regularity properties of u (see Section [3]), there is a time 
t# > such that 

(5.12) SMv\\u{t)\\HincHm<C- 

Now, let (to, ^2) C W, for some t2 > to > {k) , |to - ^21 > 1 such that (I5.12p holds 
(without loss of generality, we can assume that > t#). This claim is an immedi- 
ate consequence of the aforementioned L^-{H^ n C" (fi)) smoothing property and 
bounds ([5TTT|) . Using (|5T1|) and dEH]), we obtain 

(5.13) 



\\u {to) - u (^2)11 i2(o) = 2 I (5tu (s) , u (.s) - u (to)) ^2(0) 

|5tw(s)||(jji(t^))- (||"(s)||ffi(s:^) + ||w(io)|lHi(o)) 



/to 

< 2 



< C ||5tli||il(4„^t^.(jjl).) (||w|lL~(t.,oo;Hl) + 1) ^ 

Therefore we can choose a time t, (fc) = t < to < ^2, such that 

(5.14) h(io)-«(i2)|li2(o) < 3 

provided that (|5.7p holds for all t G (to, 12). Since G a; [uo], a large (redefined) r 
can be chosen such that 

(5.15) II"(t') "W*IIl2(o) < 3. 
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whence, (|5.14p yields (r, oo) C W. Indeed, taking 

t = inf|t >r : \\u (t) ~ u4^2(^n) > '^j , 

we have t > t and ||u (t) — w*||^2(f2) > £ if t is finite. On the other hand, in view 
of (I5.14P and (|5.15p . we have 

2(5 

||W W - U*|li2(s:j) < ||w(t)-u(T)|1^2(s:j) + |lu(T)-W,|1^2(o) < y, 

for all t > i > T, and this leads to a contradiction. Therefore, t = oo and by (|5.1ip 
the integrability of dtu in L^{t, oo; {H^ (fi)) ) follows. Hence, w [uq] = {i**} and 
(jS.ip holds on account of the L?-[H^ n C" (fi)) smoothing property. The proof is 
finished. □ 

Proof of i5.2\) . Without loss of generality, suppose now that, for all i > > 0, we 
have {u (t)) > £^ [u^,) (otherwise, there is nothing to prove). Define the function 

S(0 ■.^Ee{u{t))-Ee{u,) 

and observe that by (|5.5p - (|5.6p and (|4.23p . it satisfies 

{t) + CS {tf^^-"'^ < 0, for all t > U 
for some positive constant C = C (u) . Integration of the preceding inequality yields 

(5.16) E{t)<E (0) (l + CE (0)'-'' , 

for all t> t^,. On the other hand, we have on account (15.51) that 



{tf = -OE itf-' |s it) > ce II Vm (t)ii^.(^) , 

for all t > i*, provided that ||u (t) — u*|| < 6. Integrating this inequality over (i, oo) , 
we also get 

/oo 
l!VAi(s)||^2(o) ds<C <oo, 

for all t > t^,. As above, we obtain 



l|c'tu(s)||(^i(t^)). ds < C <oo, 
and combining with (|5.5p and (|5.16l) yields 

/oo poo 
\\^tu{s)\\^^^^^^yy ds<C \\V^l{s)\\^2(^)ds 

< C(l+t)"T^ , 

for some positive constant C, which depends on S (0) , 9, u and e. By using stan- 
dard interpolation inequalities (see. Section |4]) one can deduce the convergence 
rate estimate in the stronger norm in ()5.2p . Of course, the convergence exponent 
deteriorates. The proof is complete. □ 
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